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THE RECENT AND CONTINUING ASSEMBLY OF FIELD ELLIPTICALS BY RED MERGERS 
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ABSTRACT 

We present a study of tidal debris associated with 126 nearby red galaxies, selected from the 1.2 degree 2 
Multiwavelength Survey by Yale-Chile (MUSYC) and the 9.3 degree 2 NOAO Deep Wide-Field Survey. In the 
full sample 67 galaxies (53 %) show morphological signatures of tidal interactions, consisting of broad fans of 
stars, tails, and other asymmetries at very faint surface brightness levels. When restricting the sample to the 86 
bulge-dominated early-type galaxies the fraction of tidally disturbed galaxies rises to 71 %, which implies that for 
every "normal" undisturbed elliptical there are two which show clear signs of interactions. The tidal features are 
red and smooth, and often extend over > 50kpc. Of the tidally distorted galaxies about 2/3 are remnants and 
1/3 are interacting with a companion galaxy. The companions are usually bright red galaxies as well: the median 
/?-band luminosity ratio of the tidal pairs is 0.31, and the median color difference after correcting for the slope 
of the color-magnitude relation is —0.02 in B — R. If the ongoing mergers are representative for the progenitors 
of the remnants ~ 35 % of bulge-dominated galaxies experienced a merger with mass ratio > 1 : 4 in the recent 
past. With further assumptions it is estimated that the present-day mass accretion rate of galaxies on the red 
sequence AM/M — 0.09 ± 0.04 Gyr -1 . For a constant or increasing mass accretion rate with redshift, we find 
that red mergers may lead to an evolution of a factor of > 2 in the stellar mass density in luminous red galaxies 
over the redshift range < z < 1 , consistent with recent studies of the evolution of the luminosity density. We 
conclude that most of today's field elliptical galaxies were assembled at low redshift through mergers of gas-poor, 
bulge-dominated systems. These "dry" mergers are consistent with the high central densities of ellipticals, their 
old stellar populations, and the strong correlations of their properties. It will be interesting to determine whether 
this mode of merging only plays an important role at low redshift or is relevant for galaxies at any redshift if they 
exceed a critical mass scale. 

Subject headings: galaxies: evolution — galaxies: formation — galaxies: elliptical and lenticular, cD 



1. INTRODUCTION 

Although greatly outnumbered by spiral galaxies, elliptical 
galaxies contain ~ 20 % of the stellar mass in the present-day 
Universe (Fukugita, Hogan, & Peebles 1998) and dominate the 
high mass end of the galaxy mass function (e.g., Nakamura 
et al. 2003). Their status at the top of the food chain is a strong 
motivation for their study as it provides insight in the processes 
governing the build-up of galaxies over cosmic time, the nature 
of the most luminous galaxies at z ~ 3 and beyond, the for- 
mation of the first stars, and the origins of supermassive black 
holes. 

Despite vigorous research efforts it is still not known how 
and when elliptical galaxies were formed. Separating the time 
of the formation of their stars from the time of their assembly, 
a partial answer is that most of the stars in the most massive el- 
lipticals were formed at high redshift. This result is supported 
by a large number of studies of nearby and distant galaxies, and 
applies to ellipticals in the general field as well as those in rich 
clusters (e.g., Bower, Lucey, & Ellis 1992; Ellis et al. 1997; 
Bernardi et al. 1998, 2003; van Dokkum et al. 1998a, 1998b; 
Treu et al. 1999, 2002; van Dokkum et al. 2001a; Holden et al. 
2005). Even when accounting for possible selection effects due 
to morphological evolution ("progenitor bias"; van Dokkum & 
Franx 2001), it is very difficult to fit the observed evolution of 
massive ellipticals with formation redshifts substantially lower 
than z ~ 2. Although there is general agreement that the mean 
age is high, there are strong indications that not all stars in all 
elliptical galaxies formed at high redshift: Trager et al. (2000) 
and others (e.g., Yi et al. 2005) find evidence for a "sprin- 
kling" of recent star formation in otherwise old ellipticals, and 



recent studies of the evolution of the Fundamental Plane rela- 
tion (Djorgovski & Davis 1987) find evidence that low mass 
early-type galaxies have lower luminosity-weighted ages than 
high mass ones (Treu et al. 2005; van der Wei et al. 2005). 

Much less is known about the assembly of elliptical galax- 
ies, that is, when and how the galaxies were "put together". 
In the literature, a distinction is often made between models in 
which elliptical galaxies form from the collapse of primordial 
gas clouds ("monolithic collapse"; e.g., Eggen, Lynden-Bell, & 
Sandage 1962; Jimenez et al. 1999), or in mergers of smaller 
galaxies ("hierarchically"; e.g. Toomre & Toomre 1972; White 
& Frenk 1991). Within the framework of a merger model the 
question of assembly can be specified as, first, the typical age of 
the last major merger; and second, the nature of the progenitor 
galaxies. In many individual cases we can readily answer both 
questions: abundant evidence leaves little doubt that some el- 
lipticals formed recently through the merger of spiral galaxies, 
or will do so in the near future (e.g., Toomre & Toomre 1972; 
Schweizer 1982; Hibbard & van Gorkom 1996). However, it is 
not clear whether such mergers are exceptional or responsible 
for the formation of the majority of ellipticals. 

Among the observational arguments for recent mergers 
are the detection of "fine-structure" (e.g., shells and rip- 
ples) in a large number of nearby ellipticals (Malin & Carter 
1983; Schweizer et al. 1990); the presence of kinematically- 
decoupled cores (Franx & Illingworth 1988; Bender 1988); 
small-scale dust seen in the majority of ellipticals (e.g., van 
Dokkum & Franx 1995); (some) studies of the redshift evolu- 
tion of galaxy morphologies, and of close pairs (e.g., Le Fevre 
et al. 2000; Patton et al. 2002; Conselice et al. 2003); and 
strong evolution of the mass density of red galaxies to z ~ 1, in- 
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ferred from the COMBO-17 survey (Bell et al. 2004). Further- 
more, and perhaps most importantly, the hierarchical assem- 
bly of galaxies is a central aspect of current models for galaxy 
formation in ACDM cosmologies (e.g., Kauffmann, White, & 
Guiderdoni 1993; Cole et al. 2000; Meza et al. 2003). 

Among arguments against the merger hypothesis are the low 
probability of mergers in virialized clusters, which are domi- 
nated by ellipticals (Ostriker 1980; Makino & Hut 1997); the 
fact that the central densities of ellipticals are too high to be the 
result of dissipationless mergers of pure disks (Carlberg 1986; 
Hernquist 1992); the low scatter in the color-magnitude relation 
and Fundamental Plane, which is inconsistent with the diver- 
sity in stellar populations expected from spiral galaxy mergers 
(e.g., Bower, Kodama, & Terlevich 1998, Peebles 2002); the 
much higher specific frequency of globular clusters for ellipti- 
cals than for spirals (although some globulars may be formed 
during mergers; e.g., Schweizer et al. 1996); the observation 
that known remnants of spiral galaxy mergers have underlumi- 
nous X-ray halos (Sansom, Hibbard, & Schweizer 2000); the 
existence of the M. — a relation (Gebhardt et al. 2000; Fer- 
rarese & Merritt 2000); and the presence of massive galaxies 
at high redshift (e.g., Franx et al. 2003; Daddi et al. 2004; 
Glazebrook et al. 2004) 

These apparently contradictory lines of evidence may be 
largely reconciled by postulating that most elliptical galaxies 
formed through (nearly) dissipationless (or "dry") mergers of 
red, bulge-dominated galaxies rather than mergers of spiral 
disks. Motivated by the comparatively red colors of ellipti- 
cals exhibiting fine-structure Schweizer & Seitzer (1992) dis- 
cuss this possibility, but argue that mergers between early-type 
galaxies are statistically unlikely as the "median" field galaxy 
is an Sb spiral (see also Silva & Bothun 1998). However, most 
mergers likely occur in groups, where the early-type galaxy 
fraction is much higher than in the general field (e.g., Zablud- 
off & Mulchaey 1998). Furthermore, semianalyical models of 
galaxy formation have predicted that the most recent mergers 
of bright ellipticals were between gas-poor, bulge-dominated 
galaxies (Kauffmann & Haehnelt 2000; Khochfar & Burkert 



2003), and several studies have shown that such merging is 
consistent with the observed Fundamental Plane relation (e.g., 
Gonzalez-Garcfa & van Albada 2003; Boylan-Kolchin, Ma, & 
Quataert 2005). Finally, mergers between red galaxies may be 
common in young, unvirialized galaxy clusters at z ~ 1 (van 
Dokkum et al. 1999, 2001b, Tran et al. 2005). 

Although some early-type/early-type mergers are known to 
exist in the local universe (e.g., Davoust & Prugniel 1988; 
Combes et al. 1995), little is known about their frequency. Here 
we investigate the relevance of dry merging by analyzing the 
frequency and nature of tidal distortions associated with a well- 
defined sample of bright red galaxies. This study was motivated 
by the results of Kauffmann & Haehnelt (2000), Khochfar & 
Burkert (2003), and Bell et al. (2004), the observations of red 
mergers in z ~ 1 clusters, and the advent of very deep photomet- 
ric surveys over wide areas. We assume Cl m — 0.3, f^A = 0.7, 
and tfo = 70kms-' Mpc" 1 (Spergel et al. 2003). 

2. MORPHOLOGICAL SIGNATURES 

Mergers between gas-poor, bulge-dominated galaxies are 
generally more difficult to recognize than mergers between gas- 
rich disks. Elliptical-elliptical mergers do not develop promi- 
nent tidal tails dotted with star forming regions, but are instead 
characterized by the ejection of broad "fans" of stars, and in 
certain cases an asymmetric deformation of the inner isophotes 
(e.g., Rix & White 1989; Balcells & Quinn 1990; Combes et al. 
1995). Tails may develop if one of the progenitors rotates or 
has a disk component (Combes et al. 1995), but as there is no 
cold, young component these are expected to be more diffuse 
than those seen in encounters between late-type galaxies. 

These effects are illustrated in Fig. 1, which shows the evo- 
lution of an off-axis collision between two hot stellar systems. 
The simulation was performed with an implementation of the 
Barnes & Hut (1986) hierarchical tree code, using 65,536 parti- 
cles. The galaxies have a 3:1 mass ratio and are represented by 
Plummer models. The simulation illustrates the rapid merging 
of the two bodies, the lack of tidal tails, and the development 
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Fig. 1 . — Illustration of an off-axis collision between hot stellar systems, with mass ratio 3:1. The simulation was performed with an implementa- 
tion of the Barnes & Hut (1986) tree code, using 65,536 particles. The galaxies merge quickly; after the merger the only morphological signature 
is the presence of a broad asymmetric fan of material, which rapidly becomes more diffuse. Note the lack of tidal tails, which only form if there 
is a significant cold component. 
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of a large fan of material. The central parts of the merger rem- 
nant quickly settle in a new equilibrium configuration, whereas 
the material in the outer parts - where dynamical timescales 
are long - gradually becomes more diffuse. For a more gen- 
eral analysis of the development of elliptical-elliptical mergers, 
including the effects of rotation and different mass ratios, we 
refer to Combes et al. (1995) and subsequent papers. 

The surface brightness of tidal features associated with this 
type of mergers is comparatively low: the features consist of 
old, red stellar populations, and have higher MjL ratios than 
the blue tails associated with mergers of late-type galaxies. Fur- 
thermore, they are typically more diffuse, due to the lack of a 
cold disk in the progenitors. After the merger, the evidence 
of the past encounter becomes increasingly difficult to detect 
as dynamical evolution reduces the surface brightness further 
(e.g., Mihos 1995). No detailed simulations have been done of 
the expected surface brightness evolution of tidal features re- 
sulting from early-type mergers, but a very rough estimate can 
be obtained from the simulation shown in Fig. 1. Calibrating 
the output such that the average surface brightness within the 
effective radius of the remnant (fi R ) — 20magarcsec~ 2 (J0r- 
gensen, Franx, & Kjaergaard 1995), we find that the "fan" of 
material to the upper left has /i« <~ 25 at t = 10, i.e., immedi- 
ately after the merger. The following timesteps show a contin- 
ued rapid decrease of the brightness, reaching /j« <~ 27 at t = 15, 
the final time step of the simulation. 

It is difficult to reach such low surface brightness levels ob- 
servationally: extreme enhancements of photographic plates 
enabled Malin & Carter (1983) to detect features to /i <~ 26.5, 
similar to the limits reached by Schweizer & Seitzer (1992) 
with their deep CCD exposures on the Kitt Peak 0.9 m tele- 
scope. Detecting features with surface brightness 27 — 28 re- 
quires several hours of integration time on 4 m class telescopes, 
and to image a sizeable sample of nearby galaxies would re- 
quire hundreds of hours. Furthermore, flat fielding would have 
to be accurate to ^ 0.1 % over the entire 10' — 30' field. 

These practical difficulties can be solved by exploiting the 
fact that the (1 + z) 4 cosmological surface brightness dimming 
is a very slow function of redshift for z <C 1 . Whereas moving 
a galaxy from z = 0.01 to z = 0. 1 reduces its total brightness by 
5 magnitudes, its surface brightness changes by only 0.35 mag. 
In practice it is easier to detect very faint surface brightness fea- 
tures around a galaxy at z = 0.1 than at z = 0.01 as it occupies 
a 70 x smaller area of the detector. Furthermore, and most im- 
portantly, a study of tidal features at z ~ 0. 1 does not require a 
dedicated survey but can be done with existing imaging data. 

3. DATA 

We use data from two deep extra-galactic surveys: the Mul- 
tiwavelength Survey by Yale-Chile (MUSYC; Gawiser et al. 
2005), and the NOAO Deep Wide-Field Survey (NDWFS; 
Januzzi & Dey 1999, Januzzi et al. 2005, Dey et al. 2005). 
The data products from the optical imaging components of the 
two survey are quite similar, as they both use the 8192 x 8192 
pixel MOSAIC cameras on the NOAO 4m telescopes. The full 
dataset covers 10.5 square degrees of sky, reaching la surface 
brightness levels of [i s=a 29 mag arcsec~ 2 . 

3.1. MUSYC 

MUSYC is a deep optical/near-infrared imaging and spec- 
troscopic survey of two Southern and two equatorial 30' x 30' 

1 http://www.noao.edu/noao/noaodeep/DR3/ 



fields, chosen for their low background and the availability of 
data at other wavelengths. The fields are the Extended Chandra 
Deep Field South (E-CDFS), the Extended Hubble Deep Field 
South (E-HDFS), the field centered on the z = 6.3 Sloan Digital 
Sky Survey quasar SDSS 1030+05 (Becker et al. 2001), and the 
blank low background field CW 1255+01. The primary goal of 
the survey is the study of normal and active galaxies at z > 2. 

The survey design and data analysis procedures are described 
in Gawiser et al. (2005). Briefly, the optical imaging data con- 
sist of deep UBVRIz exposures obtained with MOSAIC-I on 
the Mayall 4m telescope on Kitt Peak and MOSAIC-II on the 
Blanco 4m telescope on Cerro Tololo. The images have a scale 
of 0."267pix~ 1 , and typical seeing of l."l. The data reach 
depths of w 26 in U, w 26.5 in B, V, and R, w 25 in /, and 
w 24 in z (AB magnitudes, 5cr point source detections), with 
some variation between the fields. The full exposure time is 
typically realized over an area of 33' x 33', and the total optical 
survey area is 1 .2 deg 2 . 

For detection of low surface brightness features we used the 
co-added "BVR" images (see Gawiser et al. 2005), which are 
combinations of the B, V, and R frames. The total effective in- 
tegration time of these frames is typically ~ 8 hours. The final 
depth reaches magnitude w 27 (AB; 5a point source detection), 
with slight variations from field to field. 

3.2. NOAO Deep Wide-Field Survey 

The NDWFS is a public optical and near-infrared imaging 
survey over two 9.3 deg 2 fields: a 3° x 3° field in Bootes, cen- 
tered at U h 32' n , +34° 17', and a 2.3° x 4° field in Cetus, cen- 
tered at 2 ;, 10 m , -4°30' (see Januzzi & Dey 1999; Januzzi et al. 
2005). The primary goal of the survey is to study large scale 
structure at z > 1 . We use the optical images available in the 
third data release (DR3; October 22 2004). This release com- 
prises the entire 9.3 deg 2 Bootes field. 

Compared to MUSYC, the NDWFS sacrifices multi-band 
coverage for area. Only three optical filters are used: B w , R, 
and /. The B w filter was designed specifically for the NDWFS: 
it is a very broad blue filter, effectively a combination of the 
standard U and B filters. Exposure times vary across the field; 
the median values are 8400 s in B w , 6000 s in R, and 1 1400 s in 
/. Details of DR3 can be found on the NDWFS webpages 1 . 

The data were not pasted together in a mosaic by NOAO but 
released as 27 partially overlapping 35' x 35' MOSAIC point- 
ings. We obtained all 27 pointings in B w , R, and / from the 
NOAO archive. In order to increase the signal-to-noise (S/N) 
ratio we created co-added "BRI" images, in the following way. 
First, all images were normalized such that 1 count corresponds 
to AB= 30 mag. Next, the summed BRI images were created 
from the B, R, and / frames according to BRI= B +R + 0.5 x I. 
This procedure optimizes the S/N ratio for objects with the ap- 
proximate SED of nearby elliptical galaxies. 

3.3. Depth of Summed Images 

The detection of diffuse tidal features depends on the limit- 
ing surface brightness level of the data, not on the point source 
detection limit. To assess the (AB) surface brightness limits 
we placed random apertures with an area of 1" in empty re- 
gions, and determined the rms fluctuations. The co-added im- 
ages were used, i.e., the "BVR" images for MUSYC and the 
"BRI" images for NDWFS. 
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Fig. 2. — Sample selection. Solid lines show the expected B — R colors and R magnitudes for L* and 3L» elliptical galaxies at redshifts 
0.04 < z < 0.20. The dashed lines shows the selection region in color and magnitude. MUSYC consists of four 0.3 deg 2 fields, including the areas 
around the Chandra Deep Field South and the Hubble Deep Field South. The NDWFS covers a contiguous area of 9.3 deg 2 . 



The MUSYC images have a typical la surface bright- 
ness limit of n w 29.5 mag arcsec~ 2 whereas this limit is w 
29 mag arcsec~ 2 for the NDWFS images. These numbers are 
approximate: the exact limit depends on the field, the local flat 
fielding accuracy, and on the assumed spectral energy distribu- 
tion of the source. The greater depth of MUSYC was expected, 
as more exposure time was devoted to the efficient B, V and R 
filters. 

The ability to detect low surface brightness emission is influ- 
enced by the extent of the features, confusion with neighboring 
objects, the contrast with the smooth emission from the galaxy, 
and other effects; in practice, we find we can confidently iden- 
tify features down to fi w 28 mag arcsec~ 2 . 

4. SAMPLE SELECTION 

The initial sample selection is based on color and magnitude 
only. No morphological criteria are applied, for three reasons: 
first, a simple color cut is straightforward to reproduce; sec- 
ond, as the galaxies are typically at z ~ 0.1 it can be difficult 
to determine morphologies from ground-based data (see § 5.1); 
and third, the classification process itself might bias the sample 
toward or against galaxies with tidal features. 

4.1. Catalogs 

The combination of the two surveys comprises 3 1 multi-band 
tiles, each covering w 0.3 deg 2 . Catalogs were created using 
the SExtractor software (Bertin & Arnouts 1996). SExtrac- 
tor catalogs were available for both surveys but these are op- 
timized for faint point sources, not for large, nearby galaxies. 
SExtractor was run with its default settings, with the follow- 
ing changes. The detection threshold was set at 20cr, with the 
added requirement that 10 adjacent pixels are above the thresh- 
old, and the mesh size for background subtraction was set at 
400 pixels rather than the default 64 to avoid oversubtraction of 



the background for large galaxies. Matched photometry in each 
band was obtained by running SExtractor in dual image mode, 
always using the /?-band images for detection. 

SExtractor's AUTO magnitudes are used as best estimates 
for total magnitudes. Colors were initially determined in fixed 
5" diameter apertures. However, we find that colors measured 
in fixed apertures are not well suited for these large objects 
as they introduce distance-dependent selection biases. Specif- 
ically, when selecting on aperture color the sample invariably 
contains large, presumably very nearby spiral galaxies with 
blue disks and red bulges. In the following, B R colors re- 
fer to Bauto — ^auto, measured in dual image mode. 

Finally, the measured B w — R colors in the NDWFS are con- 
verted to the standard Johnson system. The B w filter was cre- 
ated for the NDWFS and has a bluer central wavelength than 
the B filter. Using templates from Coleman, Wu, & Weedman 
(1980) we derive 

(B-R) = 0.27 + 0.62(B W -R)+ 0.07 (B w -R) 2 . (1) 

This transformation holds to a few percent for 1 < (B— R) < 2.5 
and redshifts z < 0.2, and was applied to the NDWFS photom- 
etry prior to the sample selection. Total magnitudes and colors 
are given on the Vega system to facilitate comparisons to previ- 
ous studies. 

4.2. Color Selection 

Galaxies are selected on the basis of their total R magnitude 
and B R color. Figure 2 shows color-magnitude diagrams for 
the four MUSYC fields (left) and for the 27 NDWFS fields 
(right). Only objects with a star/galaxy classification < 0.06 
are shown (see Bertin & Arnouts 1996). This classification 
is fairly robust at the magnitudes of interest, and effectively 
removes both non-saturated and saturated stars from the sam- 
ple. Visual inspection revealed that the software assigned a 
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star/galaxy classification > 0.06 to four galaxies in our sam- 
ple; for these objects the value was manually set to zero. We 
note that the cores of some bright galaxies can be saturated, 
which influences the magnitudes and (particularly) the colors 
in unpredictable ways. As this effect is difficult to quantify and 
only influences a handful of objects we did not attempt a cor- 
rection. The distribution of points is similar in the two surveys. 
The number of objects in the NDWFS is an order of magnitude 
larger than in MUSYC, as expected from the difference in area. 

We selected galaxies with the colors and magnitudes of 
L > L» early-type galaxies at 0.05 < z < 0.2. The two solid 
lines show the colors of the Coleman et al. (1980) E/S0 tem- 
plate redshifted from z = 0.04 to z = 0.20, normalized to 
and 3L* (Blanton et al. 2001). Dashed lines delineate the 
adopted selection region: R < 17, 1.6 < (B — R) < 2.2, and 
(B — R) > 1.6 + 0.12 x (R- 15). The galaxies in the dashed 
region are the brightest and reddest in 10.5 square degrees of 
sky. 

4.3. Additional Steps 

After visual inspection of the initial sample of 155 red objects 
32 were removed for various reasons. The discarded objects 
fall in a wide range of categories: mis-classified stars (usually 
blends with other stars or galaxies); objects on the edge of a 
field; spurious objects (always near very bright stars); severely 
saturated galaxies (although interesting in their own right, as 
they are likely very bright active nuclei); and galaxies that have 
more than one catalog listing. The latter occurs because the 27 
NDWFS tiles have some overlap; in these cases the data from 
the overlapping pointings were added, unless the S/N is much 
higher in one of the two. We also removed 12 galaxies that 
are likely members of known galaxy clusters. Eight are located 
within 8' of NSC J142841+323859, a cluster at z = 0.127, and 
four are located within 4' of the poorer, more compact cluster 
NSC J142701+341214 at z « 0.2 (Gal et al. 2003). 

Finally, we inspected all galaxies with 17 < R < 17.75 to 
identify merging pairs (see §5.2) that are split into two ob- 
jects by SExtractor, and whose combined luminosities and col- 
ors would place them in our selection region. Three such pairs 
were identified, and the brightest galaxy of each pair was added 
to the sample: 4-1190, 4-1975, and 22-2252. The final sample 
thus consists of 126 red field galaxies in MUSYC and the third 
data release of the NDWFS. 

4.4. Median Redshift 

The median R magnitude of galaxies in our sample is 16.4, 
and the median B R color is 1.92. As can be seen in Fig. 
2 these values suggest that the median galaxy is an L m 1 .2L* 
early-type galaxy at redshift z ~ 0.11. Currently we have red- 
shifts for only nine of the 123 galaxies: four in the CDF-S from 
COMBO-17 (Wolf et al. 2005), two in the SDSS 1030+05 field 
from the SDSS (York et al. 2000), and three in the CW 1255+01 
field from the SDSS. The mean redshift of these nine objects 
(z) =0.13 with a spread of 0.02, consistent with the expecta- 
tions from our color selection. 

As a further test on the selection we obtained photometric 
and spectroscopic data in ten randomly selected 3° x 3° fields 
from the SDSS. Galaxies were selected in the g — r vs. r plane 
in the same manner as in our study, using the Fukugita et al. 
(1995) transformations to convert our R and B — R limits to 
SDSS r and g — r limits. The median redshift is 0.098, with 
90% of the galaxies at z > 0.05. The rms field-to-field varia- 
tion of the median is only 0.008. 



In the following we will adopt z = 0.1 as the median red- 
shift; our conclusions change very little if we were to use, e.g., 
z = 0.08 or z = 0.13 instead. For z ~ 0.1 our observed /?-band 
limit corresponds to M R <~ — 21, and the median galaxy has 
M R ~ -22. 

5. ANALYSIS 

5.1. Morphologies 

All galaxies were assigned a morphological type by visually 
inspecting their summed images (i.e., the "BVR" images for 
MUSYC galaxies and the "BRI" images for NDWFS objects). 
The images span a large range in surface brightness levels, go- 
ing from the nearly saturated central regions to very low surface 
brightness features up to > 50 kpc away from the center. As the 
relevant dynamic brightness range is ;> 10 4 each galaxy was 
displayed at four different contrast levels simultaneously. 

The morphological types are necessarily broad. Although 
the S/N ratio is very high the spatial resolution is quite poor: 
the typical seeing is 1 ." 1 , which corresponds to 3 kpc at z = 0. 1 . 
Therefore, we only have 2 — 3 resolution elements within the 
half-light radii of many galaxies (see, e.g., J0rgensen, Franx, 
& Kjaergaard 1995). The assigned types are spiral (S), indicat- 
ing the presence of spiral arms and/or star forming regions in 
a disk; SO, indicating an early-type galaxy with an unambigu- 
ous disk component; and E/S0, indicating a bulge-dominated 
early-type galaxy. We cannot securely separate elliptical galax- 
ies from bulge-dominated or face-on SO galaxies. Contrary to 
the usual definition the E/S0 class therefore encompasses ellip- 
ticals, bulge-dominated SOs, E/SOs, and SO/Es. 

Morphological classifications are listed in Table 1, and im- 
ages of all 126 galaxies are shown in the Appendix. As ex- 
pected, the sample is dominated by early-type galaxies: of 
126 objects, 10 (8%) are classified as spirals, 30 (24%) as 
SOs (i.e., disk-dominated early-type galaxies), and 86 (68 %) 
as E/SOs (i.e., bulge-dominated early-type galaxies). Most of 
the spiral galaxies have large red bulges and faint blue arms. 
Based on their morphologies at surface brightness levels p <J 

25 mag arcsec 2 we infer that virtually all galaxies in our sample 
are red because of their evolved stellar populations, and not be- 
cause of dust (as is well known from many previous studies of 
bright red galaxies in the local Universe; see, e.g., Sandage & 
Visvanathan 1978, Strateva et al. 2001). 

5.2. Tidal Features 

Tidal features were first identified by visual inspection of the 
full sample of 126 galaxies; a quantitative analysis of distur- 
bances in the restricted sample of 86 bulge-dominated early- 
type galaxies follows in § 5.3. The flag describing tidal fea- 
tures can have one of four values: for no tidal features, 1 
for weak features, 2 for strong features, and 3 for an ongoing 
interaction with another galaxy. Galaxies in the "2" class are 
generally highly deformed merger remnants, whereas the "1" 
class indicates more subtle features. The difference is obvi- 
ously quite subjective, and in the subsequent analysis these two 
classes will generally be combined into one. A "3" classifica- 
tion implies that both the primary and secondary galaxy show 
clear distortions or tidal tails. 

The key result of our analysis is the ubiquity of tidal fea- 
tures, particularly among bulge-dominated early-type galaxies. 
In the full sample of 126 galaxies, 44 (35 %) show clear signs 
of past interactions and in an additional 23 cases (18 %) the in- 
teraction is still in progress. Only 59 galaxies (47 %) appear 
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Fig. 3. — Examples of red mergers, ordered by the progression of the interaction. The images were generated by combining the B and R frames. 
The objects are 17-596 and 17-681 (a); 19-2206 and 19-2242 (b); 1256-5723 (c); and 16-1302 (d). Panel (a) spans 5' x 5'; panels (b), (c), and (d) 
span 2.'5 x 2.'5. The tidal features are faint and red, and generally barely visible in B. Similar features are seen in a large fraction of our sample of 
123 red galaxies, in particular among the bulge-dominated early-type galaxies. Images of all objects are given in the Appendix. 



undisturbed, showing no unambiguous tidal features at the sur- 
face brightness limit of the survey. The fraction of disturbed 
objects is lowest among galaxies with a clear disk component. 
Among 40 galaxies classified as S or SO, only six (15 %) show 
evidence for past or present interactions. In contrast, among 
the 86 galaxies classified as E/SO undisturbed objects are the 
exception, as 61 (71 %) show tidal features. 

The nature and extent of the disturbances span a wide range. 
Some galaxies have clearly denned tidal tails while others show 



broad fans of stars similar to the final frames of the simulation 
shown in Fig. 1. In most cases the disturbances are subtle and 
only visible at large radii and very faint surface brightness lev- 
els, although some objects are strongly disturbed throughout. 
Often we see a mixture, e.g., a well defined tail in addition to a 
broad, smooth disturbance at very faint levels. This large range 
of properties is not surprising as it reflects the variation in the 
age of the interaction, the viewing angles, and the properties of 
the progenitors. 
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Fig. 4. — Same as Fig. 3, but now highlighting faint surface brightness levels in the summed exposures. The four systems have extensive tidal 
debris, extending to 1' or more from the center. Features can be reliably detected down to ~ 28magarcsec~ 2 . Despite their large extent the 
features contain less than 10 % of the total light of the galaxies. 



The tidal features are almost always red. In many cases the 
features are only visible in the R and / frames despite the sub- 
stantial depth of the blue exposures. They also appear smooth, 
showing no or very little evidence for clumps and condensa- 
tions. In these respects the features are very different from 
the blue, clumpy tidal tails seen in spiral-spiral interactions 
(e.g., Mirabel, Dottori, & Lutz 1992; Hunsberger, Charlton, 
& Zaritsky 1996), and from the sharp shells and ripples de- 
tected in unsharp-masked images of ellipticals (e.g., Schweizer 
& Seitzer 1992, Colbert et al. 2001). We note here that our data 



are not well suited to identify narrow features, as the FWHM 
resolution of our data is about 2 kpc at z — 0. 1 . 

The remarkable nature of these red mergers and their rem- 
nants is illustrated in Fig. 3. The figure shows two examples 
of ongoing mergers, an example of a strongly disturbed merger 
remnant, and an example of a galaxy with more subtle distor- 
tions at faint surface brightness levels, arranged in a plausible 
red merger sequence. In all cases the tidal features are smooth 
and red, and quite different from the highly structured blue tails 
associated with known mergers between gas rich disk galaxies. 
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Figure 4 shows summed images of the same objects, highlight- 
ing the faintest features and providing an indication of the sur- 
face brightness levels reached by the observations. We stress 
that these objects are fairly typical examples: as can be inferred 
from the smaller images shown in the Appendix they are by no 
means unique within our sample. 

5.3. Quantitative Identification of Tidal Features 

Quantitative characterization of the tidal features is impor- 
tant for testing the robustness of the visual classifications, mea- 
suring the flux associated with the features, assessing the ef- 
fects of changing the S/N ratio, and examining correlations be- 
tween distortions and other properties of the galaxies. Quan- 
titative criteria are also useful as a tool for future studies of 
larger samples. Standard measures of asymmetry (e.g., Abra- 
ham et al. 1996; Conselice et al. 2003) are not applicable to 
these galaxies as the features comprise only a small percentage 
of the total luminosity of the galaxies. Instead a method was de- 
veloped which determines distortions with respect to a model 
light distribution (see, e.g., Colbert et al. 2001). The method 
is only meaningful for bulge-dominated early-type galaxies, 
which make up the bulk of the sample. 

First, the galaxies are fitted by an elliptical galaxy model us- 
ing the "ellipse" task in IRAF, in three iterations. After each 



iteration a mask file is updated using the residuals from the pre- 
vious fit. The center position, ellipticity, and position angle 
are allowed to vary with radius. In cases where a second galaxy 
overlaps the primary object the two objects are fitted iteratively. 
The final model image is denoted M. Next, a "clean" galaxy im- 
age G is produced in the following way. Objects bluer or redder 
than the primary galaxy are masked, by dividing the R and B 
band images, median filtering, and identifying pixels deviating 
more than a factor two from the median color of the galaxy. To 
remove small foreground and background objects and retain the 
smooth galaxy light a "reverse" unsharp masking technique is 
used: the galaxy images are compared to Gaussian-smoothed 
versions of themselves, and pixels deviating more than a factor 
two are masked (excluding the galaxy centers). Pixels in the 
vicinity of masked pixels are also masked. Finally, a fractional 
distortion image F is created by dividing G by M. The dis- 
tortion image is convolved with a 5 x 5 median filter to reduce 
pixel-to-pixel variations. The parameter t describing the level 
of distortion is defined as 

^ | ^*x,y ^x,y \: (2) 

The tidal parameter thus measures the median absolute devia- 
tion of the (fractional) residuals from the model fit. 

The procedure is illustrated in Fig. 5, for a galaxy pair with 
no visible distortions, a galaxy with a weak tidal feature, and a 
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Fig. 5. — Illustration of the derivation of t, the quantity describing the deviations from ellipse fits to the galaxies. Three galaxies are shown: 
7-1818 (top) has no visually detected features; 18-794 (center) has a weak tidal feature, and 1256-5723 (bottom) has strong tidal features. From 
left to right are shown: the summed galaxy image; the galaxy image G with sharp features, blue features, and low S/N regions masked; the 
elliptical model fit M; and the median-filtered, noise-corrected distortion image F from which t is measured. 
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strongly disturbed object. The visually identified tidal features 
are isolated and emphasized in the distortion images F shown 
at right. The corresponding values of t vary from 8 % for the 
undistorted object to 24 % for the strongly disturbed object. 

In Fig. 6(a) the values of t are compared to the visual tidal 
classifications for all 86 bulge-dominated early-type galaxies. 
There is a strong correlation: the median value of t is 0.08 
for galaxies classified as undisturbed, 0.13 for weakly disturbed 
galaxies, and 0.19 for strongly disturbed galaxies and ongoing 
mergers. Tidal features are usually visually identified if t > 0. 1 : 
of 23 galaxies with t < 0.1 only 5 (22 %) were visually classi- 
fied as tidally distorted, compared to 56 out of 63 galaxies with 
t > 0.1 (89 %). We conclude that the visually identified distor- 
tions are robust and imply median absolute deviations from an 
ellipse fit of £ 10%. 
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Fig. 6. — Comparison of visual classifications and the tidal parameter 
t for the 86 galaxies classified as E/S0. There is a strong correlation 
between t and the visual classifications. Galaxies with visually identi- 
fied tidal features generally have t > 0.1. The histograms in (b) show 
the effect of lowering the S/N ratio by a factor ten. The majority of 
galaxies move to t < 0.1, and hence would not be identified as tidally 
distorted objects. 

The distortion maps F can be used to estimate the amount 
of star light associated with the tidal features. Visual inspec- 
tion of the distortion maps suggests that pixels deviating more 
than 15 % from the model are usually associated with visible 
features. A mask F' was created by setting pixels with values 
> 0.15 in F to 1 and all other pixels to zero. The relative flux 
in the tidal features was then estimated as follows: 

Y, F *-y*{G x , y -M x ,) 
ft — ^ 1 (3) 

with the summations over all pixels x,y. As expected, the me- 
dian value of /( is low for the 25 E/S0 galaxies with t < 0.1: 
0.01, about 1 % of the galaxy light. The median is 0.04 for 
the 63 galaxies with t > 0.1 and 0.07 for the 14 most disturbed 
galaxies with t > 0.2. From varying the cutoff in F the sys- 
tematic uncertainty in these values is estimated at <~ 30 %. We 
infer that - despite their large extent - the tidal features typi- 
cally contain only about 5 % of the total light of the galaxies. 



5.4. Ongoing Mergers 

There are 23 galaxies in the sample that show a clear tidal 
connection to a secondary object. In four cases the secondary 
object is also in the sample of 126, leaving nineteen unique sys- 
tems. The fact that tidal features have developed implies that 
dynamical friction is already at work. Assuming that the galax- 
ies were initially on nearly parabolic orbits, which is a reason- 
able assumption for field galaxies, the implication is that virtu- 
ally all these galaxies will eventually merge (J. Barnes, private 
communication). We note that our selection of merger pairs by 
tidal features is different from "standard" photometric and spec- 
troscopic selection of close pairs, where the number of mergers 
is always smaller than the number of close pairs (e.g., Patton 
et al. 2002): as expected, there are several close pairs in our 
sample which do not show tidal features and are not classified 
as ongoing mergers (see, e.g., the top panels of Fig. 5). 

The 19 merger systems are shown in Fig. 7; color images of 
these objects are in the Appendix, along with all other galaxies. 
In most cases the primary galaxy is in the sample of 126 and the 
secondary galaxy is not (because it has R > 17). In cases where 
both galaxies are in the sample, and in cases where the primary 
galaxy has a double nucleus, we designated the brightest object 
or nucleus as the primary object. The primary galaxies in the 
merging pairs are equally bright as other galaxies in the sample: 
their median R magnitude is 16.3, compared to 16.4 for the full 
sample of red galaxies. 

The merger sample comprises objects connected by a tidal 
"bridge" (e.g., 1-2874, 5-2345, 18-485/522); double nuclei in 
a highly disturbed common envelope (cdfs-374, cdfs-1100, 4- 
1 190, 1 1-962); and objects of similar brightness with disturbed 
isophotes and tidal tails or fans (e.g., 2-3070/3102, 4-1975, 
7-4247, 17-596/681, 19-2206/2242). Visual inspection of the 
merging systems indicates that the mergers are usually not be- 
tween blue, disk-dominated systems but between red, bulge- 
dominated systems. Furthermore, in many cases the secondary 
galaxy appears to be of similar brightness as the primary galaxy. 
Eight of the 126 galaxies are merging with each other: the four 
bright, red pairs are 2-3070/3102, 17-596/681, 18-485/522, and 
19-2206/2242. 

We quantified these effects in the following way. The lumi- 
nosity ratio is defined as 

h = !()(*-*)/", (4) 
L\ 

and the color difference as 

A{B-R) = (B 2 -R 2 )-{B i -R 1 ). (5) 

In the 16 cases where the primary and the secondary galaxy are 
both in the SExtractor catalog we calculated the luminosity ra- 
tio and the color difference directly (see §4.1). In three cases 
the pair is too close to be separated into two objects by SEx- 
tractor (cdfs-374, cdfs-1100, and 11-962). These objects are 
perhaps better described as single galaxies with double nuclei. 
Photometry for the nuclei was obtained from aperture photom- 
etry, with the radius of the aperture equal to half the distance 
between the nuclei. We note that the double nucleus of 11- 
962 is elongated but unresolved in the NDWFS images: the 
enlarged image shown in Fig. 7 was obtained with the OPTIC 
camera (Tonry et al. 2005) on the WIYN telescope, and has a 
resolution of Of'45 FWHM. 

As a check on the robustness of the results we also obtained 
aperture photometry in fixed 5" apertures for the 16 pairs that 
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Fig. 7. — The nineteen ongoing mergers in the red galaxy sample. The center of 11-962 is unresolved at the resolution of the NDWFS, and 
the image showing the two nuclei was obtained with the WIYN telescope in 0"45 seeing. Note that in four cases both the primary galaxy and 
its companion are in the R < 17 red galaxy sample. In three cases (4-1190, 4-1975, and 22-2252) the individual paired galaxies are fainter than 
R— 17 but their combined luminosity exceeds this limit. 
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Fig. 8. — Luminosity ratios and color differences of the 19 interacting pairs. The median luminosity ratio is 0.31, and the median color difference 
is —0.07. Panel (a) shows that the color difference correlates with the luminosity ratio, in the sense that faint companions are bluer. The broken 
line is the correlation expected for early-type galaxies on the color-magnitude relation. Color differences are plotted in (b), after removing this 
relation. The median residual color difference is —0.02, with a scatter of 0.13 magnitudes. 



are well separated. This method underestimates the luminosity 
differences, as the secondary galaxies are usually more compact 
than the primary galaxies. Nevertheless, the results are very 
similar, giving a median luminosity ratio that is only ~ 20 % 
higher than derived from the SExtractor AUTO magnitudes. 

The pair photometry is given in Table 1 . Luminosity ratios 
range from 0.04 to 0.80, i.e., from relatively minor 1:25 ac- 
cretion events to nearly equal mass mergers. The median lu- 
minosity ratio is 0.31, or a 1:3 merger. We conclude that ap- 
proximately half the ongoing interactions are major mergers. 
The median color difference is only —0.07, i.e., companions are 
typically 0.07 magnitudes bluer than the primary galaxy. As the 
companions are (by definition) fainter than the primary galax- 
ies, this small difference may be an effect of the existence of 
the color-magnitude relation. Figure 8(a) shows the relation be- 
tween color difference and luminosity ratio. There is a correla- 
tion, with the faintest companions being the bluest. The dashed 
line shows the expected correlation for elliptical galaxies on the 
(galaxy cluster) color-magnitude relation (Lopez-Cruz, Bark- 
house, & Yee 2004). In Fig. 8(b) we show the distribution of 
residual color differences after subtracting this relation. The 
median of the distribution is —0.02 and the rms scatter is only 
0.13 magnitudes. 

These results show that the ongoing mergers in our sample 
occur "within" the red sequence with very small scatter. This 
result cannot be attributed to our selection criteria: as 16 of 
the 19 pairs are separated by SExtractor there is no obvious 
bias against detecting blue companions to the primary galax- 
ies. There is a hint that the observed relation in Fig. 8(a) is 
steeper than expected from the color-magnitude relation, but 
this is largely due to the two pairs with LijL\ w 0.05 and the 
spiral/SO merger 18-485/522. Among pairs with L 2 /L x > 0.1 
whose primary galaxy is classified as early-type the residual 
scatter in A(B — R) is only 0.08 magnitudes. 



TABLE 1 

Photometry of Merger Pairs 



Object 


Li/U 


A{B-R) 


cdfs-374 


0.39 


-0.08 


cdfs-1100 


0.35 


-0.12 


cdfs-6976 


0.23 


-0.14 


1-2874 


0.05 


-0.50 


2-3070 


0.63 


-0.05 


4-1190 


0.39 


0.02 


4-1975 


0.51 


-0.16 


4-2713 


0.13 


-0.29 


5-2345 


0.04 


-0.39 


7-4247 


0.21 


-0.03 


11-962 


0.69 


-0.00 


11-1278 


0.15 


-0.07 


11-1732 


0.20 


-0.29 


14-1401 


0.24 


-0.03 


17-596 


0.65 


0.02 


18-485 


0.80 


0.22 


19-2206 


0.41 


0.08 


22-2252 


0.31 


-0.07 


26-2558 


0.20 


-0.10 



6. DISCUSSION 

6. 1 . Why Were They Missed? 

The key results of our analysis are the large number of on- 
going red mergers and the ubiquity of tidal features associated 
with (in particular) bulge-dominated early-type galaxies (galax- 
ies classified as "E/S0"). As discussed in §5.2 the large scale, 
low surface brightness features that we see are different from 
the ripples and shells that have been reported previously in 
unsharp-masked images of ellipticals at z ~ 0. Our sample is 
essentially local, in the sense that the merger rate is not ex- 
pected to evolve significantly from z ~ 0.1 to z ~ 0. Nearby 
galaxies have been studied in great detail over the past decades, 
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and before discussing the consequences of our findings we ad- 
dress the question why this preponderance of smooth red tidal 
features has not been seen before. 

The subjective nature of the visual classifications can be 
ruled out as a cause given the consistency of visual and quanti- 
tative classifications: 73 % of E/SO galaxies have t > 0.1. Small 
number statistics also play a minor role given the large number 
of objects in the sample. Given the fact that more than 90 % 
of the sample comes from the NDWFS we cannot rule out that 
the fraction of interacting galaxies is unusually high in that re- 
gion of the sky. Random 3° x 3° SDSS fields (see §4.4) show 
substantial peaks in the redshift distribution within each field, 
reflecting the fact that most red galaxies live in groups and fil- 
aments. We note, however, that the red galaxies are distributed 
rather uniformly over the NDWFS area, and that some of the 
most spectacular mergers are in the unrelated MUSYC fields. 

The most likely reason why the ubiquity of red mergers in 
the local universe was missed so far is the depth and uniformity 
of the available imaging. Classic imaging studies of nearby 
elliptical galaxies (e.g., Franx, Illingworth, & Heckman 1989; 
Peletier et al. 1990) used exposure times in the red of only 
180s- 600 s on 1 m - 2 m class telescopes, and such short ex- 
posures have been the norm ever since (e.g., Zepf, Whitmore, 
& Levison 1991; Pildis, Bregman, & Schombert 1995; Jansen 
et al. 2000; Colbert et al. 2001). Similarly, the effective expo- 
sure time of the Sloan Digital Sky Survey is about 51s per filter 
(York et al. 2000). 

Among the deepest imaging surveys are those of Malin & 
Carter (1983), who performed extreme enhancements of pho- 
tographic plates to bring out low surface brightness features, 
and Schweizer & Seitzer (1992), who used 2400 s - 3600 s in- 
tegrations with the 0.9 m telescope on Kitt Peak. These studies 
reach AB surface brightness limits of [i ~ 26.5 mag arcsec~ 2 . 
Although both studies find sharp distortions in the form of rip- 
ples and shells in a large fraction of ellipticals, their depth and 
field-of-view are not sufficient to find the large low surface 
brightness features that we report here. The exposure time that 
went into each of the images discussed here is w 27,000 s on 
4 m class telescopes equipped with modern CCDs (equivalent 
to 120 hours on a 1 m class telescope), and to our knowledge 
such long exposures have never been obtained of a significant 
sample of nearby red galaxies. 

Figure 9 shows what Fig. 4 would look like if we had exposed 
for 600 s on a 1 m class telescope. Very little remains of the dra- 
matic tidal features evident in Fig. 4, illustrating the extreme 
depth of the NDWFS and MUSYC surveys. We quantified the 
effect of the S/N ratio on the detectability of tidal features by 
artificially increasing the noise in our images by a factor of ten, 
corresponding to a decrease in exposure time of a factor of 100. 
The degraded images are of similar depth as those of Malin 
& Carter (1983) and Schweizer & Seitzer (1992), and are still 
substantially deeper than virtually all other imaging studies of 
nearby elliptical galaxies. As shown in Fig. 6(b) the distribu- 
tion of t changes substantially when the S/N is decreased. Only 
29 % of E/SOs in the degraded images have t > 0.1, compared 
to 73 % in the original images. 

We conclude that the high incidence of tidal features in our 
sample is a direct consequence of the faint surface brightness 
levels reached by the data. We also note that fiat fielding un- 
certainties may prohibit the detection of broad tidal features 
around local galaxies irrespective of exposure time. The faint 
debris in Fig. 4 has a very large extent in comparison to the 



high surface brightness regions visible in Fig. 9, and moving 
the galaxies to z ~ 0.01 would change the surface brightness 
very little but increase the sizes of the debris fields to <~ 10' or 
more. Finally, galaxy surveys in blue filters (e.g., Arp 1966) 
would classify the majority of objects in our sample as undis- 
turbed even if they met the surface brightness and flat fielding 
requirements. 
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Fig. 9. — Simulated appearance of the galaxies in Figs. 3 and 4 if 
observed for 600 s on a 1 m class telescope, a typical imaging depth 
for studies of nearby bright galaxies. Owing to its low surface bright- 
ness the large scale tidal debris is all but invisible. For comparison, 
the contour outlines the pi w 27.5 mag arcsec~ 2 level in Fig. 4. 

6.2. Recent Merger History of Bulge-Dominated Galaxies 

We first determine the implications of the observed interac- 
tions only, without applying corrections for the short duration 
of the mergers. We define a sample which contains both cur- 
rent and future bulge-dominated galaxies, consisting of the 86 
bulge-dominated galaxies (classified as E/S0), minus half of 
the six E/SOs which are interacting with each other, plus the 
disk-dominated galaxies (classified as S or SO) which are in- 
volved in a major merger. Among disk-dominated galaxies the 
only major merger is the spiral/SO pair 18-485/522, whose con- 
stituent galaxies have roughly equal luminosity. The total sam- 
ple of current and future bulge-dominated galaxies is therefore 
86 — 3 + 1 = 84. Among this sample of 84 galaxies there are 
eighteen current mergers and 41 remnants, and we infer that 
70% of current and future bulge-dominated galaxies experi- 
enced a merger or accretion event in the recent past. 

The red ongoing mergers and the red tidal features associ- 
ated with many of the E/S0 galaxies very likely sample the 
same physical process at different times. On average, galax- 
ies with strong tidal features are probably observed shortly af- 
ter the merger and galaxies with weak features are observed at 
later times. Assuming that the ongoing mergers are representa- 
tive for the progenitors of all remnants we can directly infer the 
mass ratios of the progenitors of the full sample of 59 current 
and future bulge-dominated early-type galaxies. The median 
luminosity ratio of the ongoing mergers is 0.31, and the me- 
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dian color difference is negligible after correcting for the slope 
of the color-magnitude relation. Assuming that M/L oc M° 2 
(e.g., J0rgensen et al. 1996) a luminosity ratio of 0.31 implies a 
median mass ratio of 0.23, or a 1:4 merger. 

There is also indirect evidence that the progenitors of the 
remnants were typically major mergers rather than low mass 
accretion events. Simulations by Johnston, Sackett, & Bul- 
lock (2001) show that surface brightness levels 28 < yU < 33 
are typical for the debris of small satellites such as the Lo- 
cal Group dwarfs. Similarly, the average surface brightness 
of the giant stream of M31 is /iy w 30magarcsec~ 2 (Ibata 
et al. 2001) and that of the debris from the Sagittarius dwarf 
fi v ~ 31 mag arcsec~ 2 (Johnston et al. 2001). All these val- 
ues are well beyond the detection limit of our survey. Further- 
more, the fraction of tidally disturbed galaxies is much lower 
among disk-dominated galaxies than among bulge-dominated 
galaxies. Ignoring the ongoing mergers only 8 % of disk- 
dominated galaxies show tidal features compared to 62% of 
bulge-dominated galaxies. This difference is consistent with 
the idea that the events responsible for the tidal features were 
usually sufficiently strong to disturb any dominant disk compo- 
nent. Finally, the fact that the features are typically broad and 
red suggests that the progenitors were dynamically hot systems 
with old stellar populations, consistent with the properties of 
the galaxies in the ongoing mergers. We also note that the frac- 
tion of light associated with the features (<~ 5 %) is very similar 
to the fraction of light associated with the tidal debris in the 1:3 
merger simulation shown in Fig. 1. 

We conclude that approximately 35 % of bulge-dominated 
red galaxies experienced a major merger with mass ratio > 1 : 4 
in the time window probed by our observations. This result is 
direct observational confirmation of the hierarchical assembly 
of massive galaxies. 

6.3. Merger Rate and Mass Accretion Rate 

Up to this point the analysis did not require an estimate of the 
timescale of the mergers. Such estimates are obviously uncer- 
tain, but they are necessary for turning the merger fraction into a 
merger rate and a mass accretion rate. These numbers are more 
easily compared to models and other observational studies, and 
are needed for extrapolating the results to higher redshifts. 

The merger rate can be defined in a variety of ways (see, e.g., 
Patton et al. 2002). Here it is expressed as the number of rem- 
nants that are formed per Gyr within our selection area: 



(6) 



with f m the fraction of the galaxy population involved in a 
merger - with pairs counted as single objects - and t m a char- 
acteristic timescale for the mergers. We restrict the analysis to 
the sample of nineteen ongoing mergers, as simulations of the 
fading of tidal debris around the remnants of dry mergers have 
not yet been done in a systematic way. 

Following Patton et al. (2000) we assume that the timescale 
of the mergers can be approximated by the dynamical friction 
timescale, given by 

2.64 x 10 5 r 2 v c „ 

T «°= Mm A Gyr ' (7) 

where r is the physical separation of the pairs, v c is the circu- 
lar velocity, M is the mass of the lowest mass galaxy, and In A 
is the Coulomb logarithm (see Binney & Tremaine 1987; Pat- 
ton et al. 2000). The median projected separation of the paired 



galaxies is 8."6, or 16 ±3 kpc at z = 0.10 ±0.02. Assuming 
random orientations this corresponds to a median physical sep- 
aration r = 20 ± 4 kpc. To obtain an estimate of v c we assume 
that the pairs have similar line-of-sight velocity differences as 
the seven elliptical-elliptical mergers discussed in Combes et 
al. (1995). The mean velocity difference of the Combes et al. 
pairs Av = 296±91kms- 1 (with the error determined by the 
jackknife method), implying v c = \/3Av = 513 ± 158 kms -1 
for an isotropic velocity distribution. The median mass of the 
companions is (7 ± 3) x 10 10 M Q , where we used M /Lr = 4.6 
in Solar units to convert luminosity to mass (van der Marel 
1991) and assumed z = 0.10 ±0.02. Taking In A - 2 (fol- 
lowing Dubinski et al. 1999 and Patton et al. 2000) we obtain 
r fric = 0.4±0.2Gyr. With f m = 19/122 = 0.16 ±0.03 we ob- 
tain/? = 0.4 ±0.2 Gyr" 1 . 

The effect of the mergers on the mass evolution of red galax- 
ies not only depends on the merger rate but also on the mass 
change resulting from individual mergers. The mass accretion 
rate can be approximated by 



AM/M = Rx M 2 /M 1 Gyr" 



(8) 



with M2/M1 the median mass ratio of the mergers. As shown 
in § 6.2 this ratio is approximately 0.23. With R m = 0.4 we ob- 
tain AM/M = 0.09 ± 0.04 Gyr -1 , i.e., merging increases the 
masses of galaxies on the red sequence by <~ 10% every 10 9 
years. 

For comparison to other studies it is also of interest to con- 
sider the major merger fraction within a projected separation of 
20 kpc. There are seven red pairs with luminosity ratio > 0.3 
and projected separation < 11", corresponding to a fraction of 
0.06 ± 0.02. It should be stressed that this number refers to 
mergers within the red sequence, not to the merger fraction 
within the full sample of R < 17 galaxies. The colors of the 
well-separated pairs show that red galaxies "prefer" to merge 
with other red galaxies. We have not examined the prevalence 
of mergers among luminous blue galaxies, but as discussed in 
§ 1 the stellar populations of ellipticals rule out widespread 
major mergers among this population. Therefore, the major 
merger rate in the full sample of red and blue galaxies is pre- 
sumably much lower than that within the restricted sample of 
red galaxies. A very rough estimate of the merger fraction in 
the full sample is 0.06 x -/V re d/Motai ~ 0.02, in reasonable agree- 
ment with previous studies of close pairs (see, e.g., Patton et al. 
2002, Lin et al. 2004, and references therein). 

6.4. Effects on the Evolution of the Luminosity Density 

At redshifts z < 1, the observed evolution of the luminosity 
function of red galaxies reflects passive evolution of the stellar 
populations and possible changes in the underlying mass func- 
tion. As discussed in, e.g., Mcintosh et al. (2005) these changes 
can be due to mergers, galaxies entering the red sample due to 
changes in their star formation rate, or other effects. The best 
available constraints on the evolution of the luminosity func- 
tion of red galaxies were derived by Bell et al. (2004), using 
the COMBO- 17 survey. They find that the luminosity density 
of luminous red galaxies is approximately constant out to z ~ 1, 
which is surprising given the expected evolution of a factor of 
3 — 4 in the M/L ratios of the galaxies (e.g., van Dokkum et al. 
1998a, Treu et al. 2005, van der Wei et al. 2005). A possible 
explanation is that the underlying stellar mass density evolves 
as well, compensating for passive evolution of the stellar popu- 
lations (Bell et al. 2004). 
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We first determine the effect of the observed red mergers 
only, i.e., the 52 % of red galaxies that are merger remnants 
or a merger pair. Dry mergers have no effect on the total lumi- 
nosity density, but they have a strong effect on the luminosity 
density of galaxies brighter than a fixed magnitude. The effect 
of a single generation of mergers can be approximated by 

with j(z) the luminosity density of luminous galaxies before the 
mergers and j(0) the luminosity density after the mergers. For 
f m = 0.52 and {Li/U) = 0.3 we find j(z) w 0.88 j(0). We use 
Monte Carlo simulations to test this approximation for a fidu- 
cial luminosity function with a = —0.6 and M* = —19.9 (Bell 
et al. 2004). The luminosity function is evolved backward in 
time by breaking 50 % of the galaxies into pieces with luminos- 
ity ratio 0.3. Calculating the luminosity density for M < — 19 
gives j(z) = 0.89j(0), in very good agreement with the sim- 
ple estimate given above. The conclusion is that the effect of 
the observed mergers and remnants on the luminosity density 
of bright galaxies is small, of order 10 %. 

The observed interactions only probe a relatively short pe- 
riod of at most a few Gyr. In order to extrapolate the effect 
of the mergers back in time we assume the following: 1) the 
mass accretion rate at z = 0.1 is AM/M = 0.09 ± 0.04 Gyr -1 
(see § 6.3); 2) the change in luminosity density due to mergers 
is proportional to the accreted luminosity; 3) all red galaxies 
are equally likely to undergo mergers; and 4) the mass accre- 
tion rate evolves as (1 +z) m . The value of m is treated as a 
free parameter: observational constraints on the evolution of 
the merger rate may not be quite consistent (see, e.g., Patton et 
al. 2002, Concelice et al. 2003, Lin et al. 2004), and no studies 
have specifically considered the evolution of the pair fraction 
among galaxies on the red sequence. 

The dashed lines in Fig. 10 show the predicted merger-driven 
evolution of the luminosity density of bright red galaxies with 
these assumptions, for three values of m. The model with 
m = has a constant accretion rate, and in the model with 
m =1.5 the accretion rate is 3 x higher at z = 1 than it is at 



Z = 0. The dotted lines show the evolution of the M/L B ratio 
of field early-type galaxies, as measured by van der Wei et al. 
(2005). These authors find AlnM/L B = (-1.20 ± 0.18)z for 
galaxies with M > 2 x 10 n M Q (appropriate for our sample), 
which is consistent with the independent measurement by Treu 
et al. (2005). The solid line shows the predicted evolution of j B 
when both mergers and M/L evolution are taken into account. 
Grey bands indicate the combined uncertainties in A\nM/L B 
and the merger rate. The predicted evolution of j B depends 
rather strongly on the assumed evolution of the mass accretion 
rate: it is positive for a non-evolving accretion rate, constant for 
m = 1, and negative for m > 1. 

Solid points show the luminosity density in luminous red 
galaxies as measured by Bell et al. (2004) (their Fig. 5). As 
discussed extensively by these authors the data are inconsis- 
tent with passive evolution alone (dotted curves). However, as 
can be seen in Fig. 10 the data are fully consistent with models 
that include the mass accretion rate measured here. The uncer- 
tainties in the datapoints, the measured M/L evolution, and the 
merger rate are too large to distinguish between models with 
different values of m, although the data at z > 0.8 seem to favor 
models with m > 0. 

We infer that the cumulative effects of dry mergers may be 
the dominant cause of the constant luminosity density of lumi- 
nous red sequence galaxies in the COMBO- 17 survey. There 
are several caveats: both the datapoints in Fig. 10 and the 
z = 0.1 mass accretion rate have substantial errors, leading to 
a wide range of allowed models; although the most recent gen- 
eration of mergers appears to be (nearly) dissipationless, this 
may no longer hold at z <; 0.5; and the fact that disk-dominated 
red galaxies probably evolve differently from bulge-dominated 
ones is ignored. The lack of evolution in the luminosity density 
of bright red galaxies, if confirmed, is probably due to a com- 
bination of effects with dry merging a significant, but not the 
only, contributor (see also Mcintosh et al. 2005). 
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Fig. 10. — Evolution of the B-band luminosity density of red galaxies. The dotted lines show the evolution of the M/L B ratio of field early-type 
galaxies, as measured by van der Wei et al. (2005). Dashed lines show the counter-acting effect of red mergers, with different assumptions for the 
evolution of the mass accretion rate. The mass accretion rate is 9% ±4% per Gyr at z = 0.1, and evolves as (1 +z) m . The solid line shows the 
total luminosity density. Grey bands show the combined uncertainty resulting from the uncertainties in the merger rate and the evolution of the 
M/L B ratio. Data points are from the Bell et al. (2004) analysis of the COMBO-17 dataset. There is good agreement between the models and the 
data, in particular if m > 0. 
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From an analysis of tidal features associated with bright red 
galaxies we find that <~ 70 % of bulge-dominated galaxies expe- 
rienced a merger with median mass ratio 1 :4 in the recent past. 
Expressed in other ways, <~ 35 % of bulge-dominated galaxies 
experienced a major merger involving > 20 % of its final mass, 
and the current mass accretion rate of galaxies on the red se- 
quence AM/M = 0.09 ± 0.04 Gyr -1 . Assuming a constant or 
increasing mass accretion rate with redshift it is inferred that 
the stellar mass density in luminous red galaxies has increased 
by a factor of ^ 2 over the redshift range < z < 1 . 

Neither of the two standard paradigms for elliptical galaxy 
formation appears to be consistent with our results: "mono- 
lithic" assembly at high redshift or late assembly via mergers 
of gas-rich disk galaxies. Instead, elliptical galaxies appear to 
have been assembled in mergers of bulge-dominated, red galax- 
ies. This "dry" form of merging is qualitatively consistent with 
the high central densities of ellipticals (Gao et al. 2004), their 
red colors and uniform properties, the existence of the M, — a 
relation (Wyithe & Loeb 2005), and their specific frequency of 
globular clusters. 

It remains to be seen whether widespread dry mergers are 
consistent with the slope and scatter of the color-magnitude re- 
lation (see, e.g., Bower et al. 1998, Ciotti & van Albada 2001). 
The median merger in our sample makes the brightest galaxy 
more luminous by <~ 30 % and bluer by only <~ 0.02 mag in 
B — R, and its remnant will lie within <~ 0.03 mag of the color- 
magnitude relation. However, it is doubtful whether multiple 
generations of such mergers can be accommodated, unless there 
is a strong correlation between the masses of the progenitors 
(see Peebles 2002). It will also be interesting to see whether 
the evolution is consistent with the number density of massive 
galaxies at high redshift. Again, the observed mergers only 
have a <~ 10 % effect on the mass function of red galaxies, and 
the extrapolation to higher redshifts is obviously still very un- 
certain. Furthermore, the evolution of the mass function of all 
galaxies may be different from that of the subset of galaxies on 
the red sequence. 

The high merger rate confirms predictions from hierarchical 
galaxy formation models in a ACDM universe (e.g., Kauffmann 
et al. 1993; Kauffmann 1996; Cole et al. 2000; Somerville, 
Primack, & Faber 2001; Murali et al. 2002). Furthermore, 
semi-analytical models have predicted that gas-poor mergers 
between bulge-dominated systems, rather than mergers of disk 
systems, are responsible for the formation of the most mas- 
sive ellipticals (Kauffmann & Haehnelt 2000; and, in particular, 
Khochfar & Burkert 2003). Quantitative comparisons of mass 
growth are difficult as current models do not naturally produce 
red field galaxies without star formation: the colors of field el- 
lipticals in the simulations are too blue and their M/L ratios 
too low (e.g., Kauffmann 1996; van Dokkum et al. 2001a). 
Additional mechanisms, such as heating by active nuclei, ap- 
pear to be required to halt gas cooling and star formation in 
massive galaxies (e.g., Binney 2004, Somerville 2004, Dekel 
& Birnboim 2005, Keres et al. 2005). Our results imply that 
revised models which address these issues should not only re- 
produce the "red and dead" nature of ellipticals today but also 
of their immediate progenitors - which may occur naturally if 
the progenitors have masses greater than some critical mass (see 
Cooray & Milosavljevic 2005). 

The main uncertainty in the observed merger fraction is 
the possibility that the Northern NDWFS field, which con- 
tains over 90% of the sample, is special in its frequency of 



tidally disturbed objects. This seems unlikely given its area 
of <~ 400 Mpc 2 at z = 0.1, but the issue of field-to-field vari- 
ations in the merger fraction will only be settled conclusively 
when independent fields of similar size are studied in the same 
way. The main uncertainty in the merger rate and mass accre- 
tion rate is the timescale of the mergers. Although our estimates 
broadly agree with other studies (see, e.g., Lin et al. 2004, and 
references therein), this may simply reflect the fact that similar 
assumptions lead to similar results. Modeling of red, gas-poor 
mergers has not been done in a systematic way using modern 
techniques, and it will be interesting to see what the timescales 
are for the initial coalescence and the subsequent surface bright- 
ness evolution of tidal debris. Specifically, modeling of the 19 
merging systems and 44 remnants presented here would pro- 
vide much better constraints on the merger rate and mass accre- 
tion rate, particularly when more complete redshift information 
is available. 

More detailed observational studies of the mergers and their 
remnants may help answer the question why the mergers are 
red, i.e., what made the progenitors lose their gas? If active nu- 
clei prevent the cooling of gas above some critical mass at early 
times they may play the same role during mergers, and it will be 
interesting to compare the degree of nuclear activity in undis- 
turbed galaxies, ongoing mergers, and remnants. Also, sensi- 
tive diagnostics of young populations (e.g., H<5 line strengths 
and ultra-violet photometry) can provide better constraints on 
the star formation histories of the mergers and remnants. Fi- 
nally, studies with higher spatial resolution can provide infor- 
mation on the detailed isophotal shapes (boxy or disky) of the 
remnants and their progenitors, and on possible correlations be- 
tween large scale smooth distortions and the sharp ripples and 
shells that have been reported in z ~ ellipticals (see Hernquist 
& Spergel 1995). 

Our study focuses on events that we can see today, and it will 
be very interesting to push the analysis to higher redshifts. Al- 
though the most recent generation of mergers could be largely 
"dry", previous generations likely involved blue galaxies and/or 
were accompanied by strong star formation (see, e.g., Sanders 
et al. 1988). Unfortunately it will be difficult to identify the 
broad red tidal features that we see here at significantly higher 
redshift due to the (1 + z) 4 cosmological surface brightness 
dimming. The limiting depth of the MUSYC and NDWFS im- 
ages is <~ 29 mag arcsec~ 2 (lc, AB). An equivalent survey at 
z > 1 should cover an area of ;> 1 square degree and reach lev- 
els of ~ 31.5 mag arcsec~ 2 at z = 1 and <~ 33.5 mag arcsec~ 2 
at z = 2. Even when (unfavorable) -corrections are ignored 
these requirements are well beyond the capabilities of current 
ground- or space-based telescopes. A more viable technique is 
to focus on the fraction of red galaxies in pairs. Although pair 
statistics require large corrections due to the short timescale of 
the mergers, pairs are easily detectable out to high redshift with 
the Hubble Space Telescope (see van Dokkum et al. 1999). 
Based on the work presented here the merger fraction among 
galaxies on the red sequence should be (0.06 ±0.02) x (1 +z) m 
for separations < 20kpc and luminosity ratios > 0.3. 
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APPENDIX 

CATALOG AND ATLAS 

Here we provide coordinates, magnitudes, colors, morphologies, and tidal classifications for all 126 red galaxies. We also present 
images of all objects in the sample. These small images at fixed contrast level do not do the data justice, but space limitations prohibit 
an atlas on the scale of Figs. 3 and 4. 



The atlas is available at www.astro.yale.edu/dokkum/mergers/ along with a version of the paper with much 
higher quality figures. 
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